combine 

$TITLE COMBrttp 
★ 

* GAMS RESOURCE ALLOCATION / SCHEDULING MODEL (COMBINE) 
* 

* This version combines the flow and schedule models into one integrated model 

* that accounts for the availability of personnel while optimizing flow 

* -it also uses individual rates, and has constraints to ensure a minimum amount 

* is spent in a task before someone can change tasks 
$INLINECOM { } 

SETS 

P People 

/ 

$include 1 peoplein .prn 1 ; 
/ 

PT Person Type 

/ 

$include 1 ptype . prn ■ ; 
/ 

TK Task 
:fl NO_WORK 

i3 $INCLUDE 1 tasks . prn 1 ; 

3 / 

Jj CHECKTK (TK) 

{ 1 

SU $ INCLUDE 1 checktk. prn 1 ; 

fl / 

~* TP Time Period 

L / 

*± $INCLUDE' timelS.prn 1 ; 

U 1 

'•4 S Shift Identifier 

Q $include' sdef .prn 1 ; 
/ 

ALIAS (TK,KK) 

ALIAS (TP, STRT, FIN, TPP) ; 

SETS PDATA(P, PT, S) Person to Type to Shift mapping 
/ 

$ include 1 pdata . prn ? ; 
/ 

SDATA(S, STRT, FIN) Shift Data 

/ 

$ include 1 sdat a . prn 1 ; 
/ 

PPTMAP(P,PT) person to person type map 

PTPMAP(P,TP) person to time period availability map 

CHKPTP(P,TP) defined if someone can check during time period TP 

WP(S,TP,TPP) defines min TPs that P must work in TK if starts TK in TP 

WP1 (P, TP, TPP) defines TPP's that person can work 

PSMAP(P,S) person to shift map 
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TABLE TKMAP (TK, TK) 
$INCLUDE 1 flow.prn 1 ; 



percentage flows from task to task 



SCALAR 
TPHR 
CHKHR 



time periods per hour /4/ 
minimum checking hours per pharmacist /2/ 



PARAMETERS 

RATEHR (P, TK) 
DRATEHR (TK) 
LABORHR (PT) 
LABOR (P,TK) 
RATE (P/ TK) 
RXO (TK/TP) 
MAXPTP(TP) 
MAXPTK (TK, TP) 
RXCOST(TK) 
AVAIL ( P, TP, TK) 
TPPTP(TP) 
TPAVAIL (TK, TP) 
PRFCST(P,TK) 
TIMES (TK) 
TIMEE (TK) 
NAV (TK) 
CHKELG(P) 
PAVAIL(P) 
MINTP 
COUNT 



Rate (RXs) per person per hour for each Task TK 
Rate (RXs) per person per hour for each Task TK 
Cost per person of type PT per hour 

Cost per person P per time period 
(RXs) per person per time period for each Task TK 
Initial f of Rx's in each Task (INPUT TO MODEL) 
maximum total # of people allowed during time period TP 
maximum people allowed in Task during any- time period 
RX cost per unit time in task TK 
availability of person P in Task TK at time TP 
total people in time period (AFTER SOLUTION FOUND) 
total people available in task TK during time period TP 
task preference cost per person — specified 
earliest time that someone can be scheduled to task TK 
latest time that someone can be assigned to task TK 
indicates whether task TK can be assigned at all 
indicates if a pharmacist is eligible to do a checking function 
indicates if person P is available at any time during the day 
minimum TPs before a person can change tasks 
used as a counter 



COUNT 1 

BREAK (TP, S) specifies break and lunch times for shift S 
maxtp 

total people assigned to TK during TP in solution 
number of timeperiods that a person is available to 



TPEOP(TK,TP) 
TPCHK(P) 



do checkin 



MINTPCHK(P) minimum number of time periods a person must check 



*** The following figures out TP availability of shifts 
SETS 

STRTTP(TP,S) starting time period for shift S 
FINTP(TP,S) finishing time period for shift S 
STPMAP(TP,S) shift to time period map 



PARAMETERS 
START (S) 
FINISH(S) 



order (sequence) of starting time period for shift S 
order (sequence) of finishing time period for shift S 



BREAK (TP, S)=0 ; 

$ INCLUDE 1 SPARAMS . PRN 1 ; {file has lunch and break information for each shift} 
STRTTP(TP,S) = SUM (FIN, SDATA ( S , TP, FIN) ); 
FINTP(TP,S) = SUM(STRT, SDATA ( S , STRT , TP) ); 
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LOOP( (TP,S) $STRTTP(TP,S) , 
START (S) = ORD(TP) ; 

); 

LOOP ( (TP, S) $FINTP (TP/ S) , 
FINISH(S) = ORD(TP) ; 

) ; 

STPMAP (TP, S) $ (ORD (TP) GE START (S) AND ORD(TP) LE FINISH (S) 

AND NOT BREAK (TP, S) ) = 1; 
STPMAP (TP, S)$ (START (S) EQ FINISH (S) ) = 0; 
*** *** 

PTPMAP(P,TP) = SUM ((PT,S)$PDATA(P,PT,S), STPMAP (TP, S) ) ; 
PSMAP(P,S) = SUM(PT$PDATA(P,PT,S) , PDATA ( P, PT , S ) ) ; 
PRFCST(P, 1 NO_WORK f ) =1; 

$ INCLUDE 1 PREF.PRN' ; {file describes preference list for each person} 
AVAIL (P, TP, TK)$PRFCST(P,TK) = PTPMAP ( P, TP) ; 



THESE ARE INITIALIZATIONS STATEMENTS 
AND MUST BE BEFORE 

PARAMS . PRN 
IS READ 



IMES(TK)=1; 

IMEE (TK) =CARD (TP) ; 

AV(TK)=0; 
MAXTP = 96; 
MINTP = 1; 
MAXPTK (TK, TP) =900; 

$INCLUDE f PARAMS . PRN ' ; { also includes times and timee information] 
AVAIL(P,TP,TK)$(NAV(TK) OR ORD(TP) LT TIMES (TK) OR ORD(TP) GT TIMEE (TK) ) = 0 
MAXPTK (TK, TP) $ (NAV (TK) OR ORD(TP) LT TIMES (TK) OR ORD(TP) GT TIMEE (TK)) = 0; 
TPAVAIL (TK, TP) = SUM (P, AVAIL (P, TP, TK) ) ; 
MAXPTK (TK, TP) =MIN (MAXPTK (TK, TP) , TPAVAIL (TK, TP ) ) ; 

* THE FOLLOWING FIGURES OUT TP GROUPINGS FOR MINIMUM TIME IN TASK 

LOOP(S, 
COUNT = 0; 

LOOP(TP$STPMAP(TP,S) , 
I F (MOD ( COUNT , MINTP ) =0 , 
COUNT 1 = 0; 

LOOP (TPP$ (STPMAP (TPP, S) AND ORD(TPP) GT ORD (TP) ) , 
COUNT 1 = COUNT 1 + 1; 

IF((COUNTl LT MINTP), WP(S,TP,TPP) = YES ); / 
) ; 

IF((COUNTl LT (2*MINTP-1) ) , 

WP (S,TP,TPP) $ (STPMAP (TPP, S) AND ORD (TPP) GT ORD (TP) ) =YES ; 

) ; 
); 

COUNT = COUNT + 1; 

); 
); 

WPl (P, TP, TPP) = SUM(S$PSMAP(P,S), WP ( S , TP, TPP) ) ; 
display wp,wpl; 

★ 

* THE FOLLOWING COULD BE SPECIFIED IN THE INTERFACE/EXTERNAL FILE 

LABORHR (PT) = 11; {initialize labor cost} 

LABORHR ( 1 Clerical 1 ) = 7.5; 
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LABORHR( 'Technician 1 ) = 11; 
LABORHR ( 1 Pharmacist 1 ) = 27; 



PPTMAP(P,PT) - SUM(S$PDATA(P / PT / S) , PDATA ( P, PT , S ) ); 

PAVAIL(P) = SUM( (TP,TK)$AVAIL(P,TP,TK) , AVAIL ( P, TP, TK) ) ; 

LABOR (P,TK) = SUM(PT$PPTMAP(P,PT) , LABORHR ( PT ) ) / TPHR; 

LABOR (P, 'NO_WORK' ) = LABOR ( P / 1 NO_WORK ' ) - ( . 1* (1- . 01*CARD (TK) ) ); 

$ INCLUDE ' PRATE . PRN 1 ; 

RATEHR ( P, TK) $ ( PRFCST ( P, TK) AND RATEHR (P,TK) = 0) = DRATEHR (TK) ; 
RXCOST(TK) = 1-.01*ORD(TK) ; 
MAXPTP(TP) = 900; 

RATE (P, TK) = RATEHR (P,TK) / TPHR; 

CHKELG(P)$PPTMAP(P, 'Pharmacist' ) = SUM (CHECKTK(TK) , PRFCST (P, CHECKTK) ) ; 
CHKPTP(P / TP)$CHKELG(P) = SUM (CHECKTK, AVAIL (P, TP , CHECKTK) ) ; 
TPCHK(P)$CHKELG(P) = SUM (TP$CHKPTP ( P, TP) , CHKPTP (P, TP) ) ; 
^DISPLAY TPCHK; 
jMINTPCHK(P) = MIN (TPCHK (P) , TPHR * CHKHR) ; 

jPOSITIVE VARIABLES 

\ RX(TK,TP) # of Rx's in task TK during time period TP 

RXPROC(TK,TP) # of Rx's PROCESSED in task TK during time period TP 
SCHED (P, TP, TK) equals 1 if person P is assigned to task TK during time perio 
1 d TP 



* INTEGER VARIABLES 

-VARIABLE 
J / OBJ; 

|/equations 

COST objective function 

RXPRC(TK,TP) calculates # of Rxs processed during time period TP 
RXPRC1 (TK, TP) added condition to calculate # of Rxs processed during TP 
RXQ (TK, TP) calculates total # of Rxs in task TK during time period TP 
MXPTK (TK, TP) maximum # of people in task 

SUPPLY (P, TP) each person can only be assigned 1 task at a given time 
TEAMPC (TK, TP) team functions must have the same # of pharmacists as clericals 
CHKMIN(P) each pharmacist must check at least 2 hrs a day 

TKMIN(P,TP,TPP,TK) forces person in task to stay for min TPs before changing 
TKMAX (P, TK) maximum number of time periods a person can do a task 



COST.. OBJ =E= SUM( ( P, TP, TK) $AVAIL ( P, TP, TK) , 

(LABOR (P, TK) + PRFCST (P, TK) ) * SCHED ( P, TP, TK) ) 
+ SUM ( (TK, TP) $ (ORD (TK) <CARD (TK) ) , RXCOST (TK) *RX (TK, TP) ) ; 



*sched.up(p,tp,tk)$(NOT AVAIL (P, TP, TK) ) = 0; 
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*RXPROC.FX(TK,TP)$ (MAXPTK (TK, TP ) = 0) = 0; 

RXPRC(TK,TP) . . RXPROC(TK,TP) =L= SUM ( P$AVAIL ( P, TP, TK) , 

SCHED(P,TP,TK) * RATE (P, TK) ) ; 
RXPRC1(TK,TP) . . RXPROC(TK,TP) =L= RX (TK, TP) ; 

RXQ (TK, TP) . . RX(TK,TP) =E= RX(TK,TP-1) +RX0(TK,TP) + 

SUM (KK$TKMAP (KK, TK) , TKMAP (KK, TK) *RXPROC (KK, TP-1 ) ) 
- SUM (KK$ TKMAP (TK, KK) , TKMAP (TK, KK) *RXPROC (TK, TP-1 ) ); 
MXPTK (TK, TP) . . SUM(P$AVAIL(P,TP,TK) , SCHED ( P, TP, TK) ) =L= MAXPTK (TK, TP) ; 
SUPPLY (P, TP) $PTPMAP(P, TP) . . SUM (TK$AVAIL ( P, TP , TK) , SCHED ( P, TP, TK) ) =E= 1 ; 

* the following constraints are "hard coded" into the program -(5/6/97)-** 

TEAMPC (TK, TP) $ (CHECKTK (TK) AND MAXPTK (TK, TP) GT 0).. 

SUM(P$PPTMAP(P, 'Pharmacist') , SCHED (P, TP, TK) ) - 
SUM(P$PPTMAP(P, 'Clerical' ) ,SCHED(P, TP, TK) ) =E= 0; 

JHKMIN(P) $ (CHKELG(P) and PAVAIL (P) ) . . 

SUM ( (TP, CHECKTK) $ (PRFCST(P, CHECKTK) and AVAIL (P, TP, CHECKTK) ) , 
SCHED (P, TP, CHECKTK) ) =G= MINTPCHK ( P) ; 

Ik the following constraints ensure that people stay in task for minimum time 

[tKMIN (P, TP, TPP, TK)$ (PAVAIL (P) AND WP1 (P, TP, TPP) 

AND PRFCST(P,TK) AND AVAIL ( P, TP, TK) AND AVAIL ( P, TPP, TK) ). . 
SCHED (P, TP, TK) - SCHED (P, TPP, TK) =E= 0; 

TKMAX(P,TK)$ (PAVAIL (P) AND PRFCST (P, TK) ) . . 

SUM(TP$AVAIL(P,TP,TK) , SCHED ( P , TP, TK) ) =L=MAXTP; 

$OFFSYMXREF OFFSYMLIST 
MODEL PSCHED /ALL/; 
OPTION ITERLIM = 950000; 
OPTION* RESLIM = 54000; 
^OPTION LIMROW =3; 
OPTION LP = CPLEX; 
OPTION optca = 100; 
*PSCHED.optfile =1; 

SOLVE PSCHED USING LP MINIMIZING OBJ; 
TPEOP(TK,TP)=SUM(P, SCHED . L ( P , TP , TK) ); 

FILE REP1 / PHARM3A1.TXT /; 

FILE REP2 /PEOPLE. prn/; 

FILE REP3 /RXPROC.prn/; 

FILE REP4 /RXCUM.prn/; 

PUT REP1; 
REP1.PW=500; 

PUT "OBJECTIVE VALUE":18, OBJ.L:10:2 /; 

PUT "ABSOLUTE GAP ": 18 , PSCHED . OPTCA: 15 : 3 /; 

PUT "RELATIVE GAP " : 18 , PSCHED . OPTCR : 8 : 3 //; 
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PUT "TIME "," "," RX IN "," RX'S "," "/; 

PUT "PERIOD " , "TASK "," QUEUE "/'PROCESSED"," PEOPLE"/; 
LOOP( (TP, TK) $RX.L(TK,TP) ", 

PUT TP.TL:8, TK.TL:8, RX . L (TK, TP) : 8 : 1 , RXPROC . L (TK, TP) : 10 : 1 , 

TPEOP(TK,TP) :8:1 /; 

) ; 

PUT //; 

PUT "TIME ", "TOTAL " /; 
PUT "PERIOD "/'PEOPLE" /; 

LOOP (TP, 

TPPTP(TP) = SUM(TK, TPEOP(TK,TP) ) 
PUT TP.TL:8, TPPTP(TP):8 /; 

); 

PUT /; 



<PUT REP2; { PERSONNEL ALLOCATION REPORT } 

-REP2.PW-255; 

= *REP2.PC = 6; 

"^PUT "* PERSONNEL ALLOCATION PER TIME PERIOD" /; 

PUT " " ; 

: LOOP(TK$ (ord(TK)<18) , 
PUT TK.TL:11; 

); 

PUT /; 
LOOP (TP, 
PUT TP.TL:4; 

LOOP(TK$(ord(TK)<18) , 

PUT ROUND (TPEOP(TK, TP) +0.499) :11:1; 

) ; 
PUT /; 

); 

put " + " /; 

PUT " 

LOOP(TK$(ord(TK)>=18 AND ORD (TK) <CARD (TK) ) , 
PUT TK.TL: 11; 
) ; 

PUT /; 

LOOP (TP, 

PUT TP.TL:4; 

LOOP(TK$ (ord(TK)>=18 AND ORD (TK) <CARD (TK) ) , 
PUT ROUND (TPEOP(TK, TP) +0.499) :11:1; 

); 

PUT /; 

); 

PUT /; 
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PUT REP3; 
REP3.PW=500; 

PUT " # OF RX'S PROCESSED IN EACH TASK DURING EACH TIME PERIOD" /; 
PUT " 

LOOP(TK$ (ord(TK)<18) , 
PUT TK.TL:11; 

); 

PUT /; 
LOOP (TP, 
PUT TP.TL:4; 

LOOP(TK$(ord(TK)<18) , 

PUT RXPROC.L(TK,TP) :11:1; 

); 

PUT /; 
); 

PUT /; 

PUT " "; 

LOOP{TK$ (ord(TK)>-18 AND ORD (TK) <CARD (TK) ) , 
PUT TK.TLrll; 

); 

PUT /; 
LOOP (TP, 
PUT TP.TL:4; 

LOOP(TK$ (ord(TK)>=18 AND ORD (TK) <CARD (TK) ) , 

PUT RXPROC.L(TK,TP) :11:1; 

); 

PUT /; 

); 

PUT /; 

PUT REP4; 
REP4.PW=500; 

PUT "TOTAL # OF RX'S IN QUEUE FOR EACH TASK DURING EACH TIME PERIOD" / 
PUT " 

LOOP(TK$ (ord(TK)<18) , 
PUT TK.TL:11; 
) ; 

PUT /; 
LOOP (TP, 
PUT TP.TL:4; 

LOOP(TK$ (ord(TK)<18) , 

PUT RX.L(TK,TP) :11:1; 

); 

PUT /; 
) ; 

PUT /; 
PUT " 

LOOP(TK$ (ord(TK)>=18 AND ORD (TK) <CARD (TK) ) , 
PUT TK.TL:11; 
) ; 

PUT /; 
LOOP (TP, 
PUT TP.TL:4; 
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L00P(TK$(ord(TK)>=18 AND ORD (TK) <CARD (TK) ) , 

PUT RX.L(TK,TP) :11:1; 

) ; 
PUT /; 

) ; 

PUT /; 



$ INCLUDE ' PSCHEDFT . RPG 1 

* file testrep /test. out/; 
file DUMMY /! PSCHED/; 

* put testrep; 

* loop( (p,tk,tp)$sched.l(p,tp,tk) , 

* put p.tl:5, tp.tl:6, tk.tl:12, sched. 1 (p, tp, tk) /; 

* >; 



put DUMMY; 
PUT PSCHED. 
Q PUT$ (PSCHED, 
PUT$ (PSCHED. 
PUT$ (PSCHED, 
PUT$ (PSCHED, 
PUT$ (PSCHED. 
PUT$ (PSCHED, 
PUT$ (PSCHED, 
PUT$ (PSCHED, 
PUT$ (PSCHED 
PUT$ (PSCHED 
PUT$ (PSCHED 
PUT$ (PSCHED 



MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 
MODELSTAT 



/; 

EQ 1) "OPTIMAL"/; 

EQ 2) "LOCALLY OPTIMAL"/; 

EQ 3) "UNBOUNDED"/; 

EQ 4) "INFEASIBLE"/; 

EQ 5) "LOCALLY INFEASIBLE"/; 

EQ 6) "INTERMEDIATE INFEASIBLE"/; 

EQ 7) "INTERMEDIATE NONOPTIMAL" / ; 

EQ 8) "INTEGER SOLUTION"/; 

EQ 9) "INTERMEDIATE NON-INTEGER"/; 

EQ 10) "INTEGER INFEASIBLE"/; 

EQ 12) "ERROR UNKNOWN"/; 

EQ 13) "ERROR NO SOLUTION"/; 
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pschedft . rpg 

* REPORT GENERATOR FOR GANNT CHART (IN GAMS) 

FILE REP99 /ASSIGN1 . PRN/ ; 

PARAMETER COUNT; 

PUT REP99; 
REP99.PW = 500; 

PUT "OBJECTIVE VALUE":18, OBJ.L:10:2 /; 

PUT "ABSOLUTE GAP " : 18 , PSCHED . OPTCA: 15 : 3 /; 

PUT "RELATIVE GAP " : 18 , PSCHED. OPTCR: 8 : 3 //; 

PUT " PERSONNEL ALLOCATION PER TIME PERIOD" /; 
PUT " 
LOOP (TP, 

IF(ORD(TP) = 44, PUT / ) ; 

PUT TP.TL:5 

) ; 

PUT /; 

OLOOP ( (P,PT)$PPTMAP(P,PT) , 

^PUT P.TL:11,P.TE(P) :19, " ", PT.TL: 11; 

W LOOP (TP, 

« s COUNT = 0; 

=J IF(MOD(ORD(TP) , 19) = 0, PUT / ); 

^ IF (NOT PTPMAP(P,TP) , PUT "N/A " : 11; 

iH ELSE 

J LOOP (TK$SCHED.L(P,TP,TK) , 

COUNT = COUNT + 1; 
;=! IF (COUNT = 1, 

in PUT TK.TL:11; 

^ ) ; 
); 
) ; 
T ); 

H-PUT /; 

); 
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